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Abstract 

Somos 4 sequences are a family of sequences defined by a fourth- 
order quadratic recurrence relation with constant coefficients. For par- 
ticular choices of the coefficients and the four initial data, such recur- 
rences can yield sequences of integers. Fomin and Zelevinsky have 
used the theory of cluster algebras to prove that these recurrences 
also provide one of the simplest examples of the Laurent phenomenon: 
all the terms of a Somos 4 sequence are Laurent polynomials in the 
initial data. The integrality of certain Somos 4 sequences has previ- 
ously been understood in terms of the Laurent phenomenon. However, 
each of the authors of this paper has independently established the 
precise correspondence between Somos 4 sequences and sequences of 
points on elliptic curves. Here we show that these sequences satisfy a 
stronger condition than the Laurent property, and hence establish a 
broad set of sufficient conditions for integrality. As a by-product, non- 
periodic sequences provide infinitely many solutions of an associated 
quartic Diophantine equation in four variables. The analogous results 
for Somos 5 sequences are also presented, as well as various examples, 
including parameter families of Somos 4 integer sequences. 

1 Introduction 

The study of integer sequences generated by hnear recurrences has a long 
history, and, apart from their intrinsic interest in number theory, nowa- 
days there are many apphcations of such sequences in computer science and 
cryptography. However, not so much is known about sequences generated 
by nonlinear recurrences. In [7] (see e.g. section 1.1.20) it is suggested that 
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Somos sequences, which are generated by quadratic recurrence relations of 
the form 

i=i 

(where the coefficients aj are constant), are suitable generalizations of linear 
recurrence sequences, in the sense that they have many analogous properties. 
In this paper we are mainly concerned with the cases k = 4 and A; = 5, 
and the problem of determining when such recurrences yield sequences of 
integers; unlike the situation for linear recurrences, this problem is not so 
straightforward . 

Somos 4 sequences are defined by a fourth-order quadratic recurrence 
relation of the form 

An+^An = a An+3An+l + (3 (A„+2)^ (1-1) 

where a, /3 are constant coefficients. We also refer to equations like (II. ip 
as bilinear recurrences, by analogy with the Hirota bilinear form of inte- 
grable partial differential equations in soliton theory [12], which have dis- 
crete (difference equation) counterparts [36]. During an exploration of the 
combinatorial properties of elliptic theta functions, Michael Somos intro- 
duced the sequence defined by (jl.ip with coefficients a = /? = 1 and initial 
data Ai = A2 = A3 = A4^ = 1, which begins 

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, . . . (1.2) 

(and in this case the sequence extends symmetrically backwards as . . . , 59, 23, 
7, 3, 2, 1,1,1,1,2,...). It is a remarkable fact that although the recurrence 
defining the sequence (jl.2p is rational (because one must divide by An in 
order to obtain An+A), this sequence consists entirely of integers. In what 
follows we are mainly concerned with Somos 4 sequences with integer coef- 
ficients and initial data, and we establish a broad set of sufficient conditions 
for integrality of such sequences. However, it is most convenient to work 
with sequences over C and indicate when our results reduce to sequences 
over Q or to integer sequences. Given four adjacent non-zero values, the re- 
currence (jl.ip can be iterated either forwards or backwards, so it is natural 
to define the sequences for n G Z. (In fact, if a zero is encountered then the 
sequence can be analytically continued through it, either by using the fact 
that the recurrence (II. ip has the singularity confinement property of [lOj, or 
by making use of the explicit formula for the iterates given in (jl.lip below.) 

The particular sequence (jl.2p above is sometimes referred to as the Somos 
4 sequence ^29j , and the first published proof of its integrality appears in the 
survey by Gale [S]. However, following [53] we will refer to (|1.2p as the Somos 
(4) sequence, in order to distinguish it from other Somos 4 sequences defined 
by recurrences of the general form (jl.ip . Such sequences originally arose in 
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the theory of elliptic divisibility sequences (EDS), which were introduced 
by Morgan Ward [331 El]) who considered integer sequences satisfying a 
recurrence of the form 

Wn+AWn = {W2f Wn+Z Wn+l " Wi (t^n+s)', (1-3) 

these being a special case of (jl.ip with 

a = {W2f, P = -WiW3. (1.4) 



The form of the recurrence (II. 3p defines an antisymmetric sequence, so that 
Wri = —W-r,. Ward showed [Ml that if the initial data are taken to be 



Wo = 0, Wi = l, W2,W3,W4eZ with VF2IW4, (1.5) 

then it turns out that the subsequent terms of the sequence are all integers 
satisfying the divisibility property 

Wn\Wm whenever n\m. (1-6) 

In this sense, EDS constitute a nonlinear generalization of Lucas sequences. 
Furthermore, Ward demonstrated that the terms of a generic EDS defined 
by ()1.3p correspond to the multiples [n]P of a point P on an elliptic curve 
E, and the nth term in the sequence can be written explicitly in terms of 
the Weierstrass sigma function associated with the curve E, as 

cr(Kj" 

Somos 4 sequences, and EDS in particular, are of considerable interest 
to number theorists due to the way that primitive prime divisors appear 
OE]- In that and in other respects, they have a lot in common with linear 
recurrence sequences (see [7], chapters 10 & 11). Somos 4 sequences have 
also arisen quite recently in a different guise, providing one of the simplest 
examples of the combinatorial structures appearing in Fomin and Zelevin- 
sky's theory of cluster algebras. More precisely, the following result was 
proved in [8j. 

Theorem 1.1. The Laurent property for Somos 4 sequences: For a 

Somos 4 sequence defined by a fourth- order recurrence U.l\) with coefficients 
a, (3 and initial data ^1,^2,^3,^4, all of the terms in the sequence are 
Laurent polynomials in these initial data whose coefficients are in 'L[a,f5], 
so that An G Z[a,/3, Af^] for all n e Z. 

The above result is only one of a series of analogous statements for a wide 
variety of recurrences in one and more dimensions to which the machinery of 
cluster algebras was applied in |8| . The integrality of the Somos (4) sequence 
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()1.2p . and also of any other sequence defined by the recurrence (ll.ip with 
coefficients a, /3 G Z and initial values = yl2 = ^3 = ^4 = 1, is an 
immediate consequence of the Laurent property. However, there are many 
other examples of Somos 4 sequences that consist entirely of integers, but 
for which the Laurent phenomenon alone is insufficient to deduce integrality. 
For instance, consider the Somos 4 sequence generated by 

An+iA^ = 1331A„+3^n+i + 119790^2^2 (1-8) 

with the initial values 

Ai = 1, A2= 3, A3 = 121, Ai = 177023. 

This turns out to be an integer sequence, i.e. An E Z for all n E Z, but 
Theorem 11.11 is not enough to show this; it follows from Theorem 13.11 be- 
low, which relies on the connection between Somos 4 sequences and elliptic 
curves. 

In the rest of the paper we shall make use of the fact that the terms A^ 
in a Somos 4 sequence correspond to a sequence of points Q + [n]P on an 
associated elliptic curve which was proved by each of us independently, 
in [21] and |13| respectively. Our results can be combined and summarized 
as follows. 

Theorem 1.2. If {An) is a Somos 4 sequence defined by a recurrence 
with coefficients a, (3, then the quantity 



« A — + ^r-i — ] + -A -A — (1-9) 



AnAn+l \An-lAn+l An An+2 J ^n-l An+2 

is independent of n, and for a 7^ the sequence corresponds to a sequence 
of points Q + [n]P on an elliptic curve E with j-invariant 

(T4-8/3T2-24a2r + 16/32)3 



' a^{f)T^ + a2r3 - 8/32r2 - 36a2/?r + 16/33 - 27a4) 
Moreover, the general term of the sequence can he written in the form 

A. = ai"'J2l±^, (1.11) 

where a{z) = (t(z; (72,53) denotes the Weierstrass sigma function associated 
with the curve E written in the canonical form 

E: y2 ^ 4rc3 - g^x - 53, (1-12) 

and the six parameters g2, g3,a,b, zq, k E C are determined uniquely (up 
to a choice of sign) from the two coefficients a, (3 and the four initial data 
Ai, ^2, A3, 744. The (j;, y) coordinates of the corresponding sequence of points 
on the curve E are given by Q + [n]P = ( p{zQ + hk) , p'{zo + hk) ). 
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Remarks. Given four non-zero initial data ^1,^2)^3)^4; the quantity T 
in (jl.9p can be calculated as 



^1712^3^4 

(For a discussion of zeros in Somos sequences, see [31].) The translation 
invariant was first found in [3T], and also appeared in [H], where it was 
denoted J, and written in terms of fn = An-i An+i/{An)'^ as 

T = fnfn+i + a(^ + -^)+-^. (1.14) 



fn fn+1 J fn fn+1 

(Note that the equation T = constant defines a quartic curve in the (/„, fn+i) 
plane: this curve is birationally equivalent to the curve E given in (11.120 .) In 
the paper [13], a more complicated invariant was used, namely the quantity 

A = (rV4-/3)/(3a). (1.15) 

We can also use the equation (jl.ip to write T as a more compact expression 
in terms of five adjacent terms of the sequence, that is 

An-l An+2 ^ aAl An-2 ,^ ^ 

= ~~x — 1 1" ~x 1 1 — 1 1 — (l-lo) 

■^n -^n+l -^n—1 -^n+1 ^n—1 -^n 

The explicit formulae for the invariants §2, 53 of the curve ()1.12p in terms 
of a, (3 and T are given in section 2 below, along with a summary of other 
results that are needed here. It is worth mentioning that the above theorem 
still holds in the degenerate case when the bracketed expression in the nu- 
merator of (jl.lOp vanishes, so that j = 00 and the curve E becomes singular; 
in that case the formula (jl.lip is written in terms of hyperbolic functions 
(or, when g2 = 93 = 0, in terms of rational functions). 

The connection between Somos 4 and Somos 5 sequences and elliptic 
curves was previously understood in unpublished work of several number 
theorists: see Zagier's discussion of Somos 5 [37], the results of Elkies quoted 
in [2], and the unpublished work of Nelson Stephens mentioned in |31| . 
where alternative Weierstrass models for the elliptic curve E are used. For 
another approach using birationally equivalent quartic curves and continued 
fractions, see |19j . A brief history of Somos sequences can be found on 
Propp's Somos sequence site [22] . 



2 Companion elliptic divisibility sequences 

The purpose of this section is to present some facts about Somos 4 sequences 
and their relationship with elliptic curves and division polynomials which 
will be useful for what follows. Given a Somos 4 sequence with coefficients a, 
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/? and initial data Ai, A2, A^, A^we can immediately compute the translation 
invariant T given by (jl.9p . and then the invariants for the associated curve 
()1.12p can be calculated directly from the formulae 

- 8/3r2 - 24a2r + 16/?2 
9^ = ^2 ' (2-1) 

- 12/5r4 - 36a2r3 + A8(3'^T'^ + UAa'^(3T + 2160^ - 64/33 
= 216^3 ' (2-2) 

while j = 172852 / idi ~ 27(?3) yields the expression ()1.10p for the j'-invariant. 
(Upon substituting for A as in p.lSp . the formulae (|2.ip and (|2.2p are seen 
to be equivalent to the formulae given in |13j.) 

The backward iterates Aq, A^i and A^2 are obtained by applying (jl.ip 
in the reverse direction, so that the quantities 

/_i = A_2 ^0/(^-1)', fo = A^iA,/iAof, /i =^0^2/(^1)', 

can be calculated, and then the associated sequence of points on the curve 
(|1.12p is given by 

Q + [n]P = (p(zo), P'(^o)) +M (p(k), p'(k)) 

(2.3) 

= (a - /o, {fo)Hfi - + N (a, ^/h) 

with A as in (jl.lSp . Up to the overall ambiguity in the sign of ^/a = p'(k), 
which corresponds to the freedom to send k —k, zq — > —zq on Jac(-E'), 
the coefficients and the translation invariant are given as elliptic functions 
of K by 



a 



p'{Kf, /? = p'(^)2(p(2k)-p(k)), T = p"{k). (2.4) 



Elliptic divisibility sequences (EDS) were originally defined by Morgan 
Ward [33] as sequences of integers satisfying the property (II. 6p as well as 
the family of recurrences 



—m 



WmWn-1 Wm-lWn 
Wm+lWn WmWn+1 



(2.5) 



for all m, n, and they were called proper EDS if Wq = 0, Wi = 1 and 
W2W3 / 0. Here we define a generalized EDS to be a sequence { Wn }nez 
specified by four non-zero initial data 

Wi = l, W2,W3,W4eC* 

together with the recurrence p.3p . of which immediate consequences are 
that Wq = and the sequence is consistently extended to negative indices 
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so that W-n = —Wn and ()1.3p holds for ah n G Z. We may then make 
the observation that a (generaUzed) EDS is just a special type of Somos 4 
sequence, and thus it follows immediately from Theorem 1.2 (i.e. the formula 
(jl.lip with a = b = 1 and zq = 0) that the terms of the EDS are given by 
the explicit expression ()1.7|) in terms of the Weierstrass sigma function. The 
whole family of recurrences ()2.5p then follows as a direct consequence of 
this explicit formula together with the three-term equation for the sigma 
function (see §20.53 in [35]). If we further assume that we have integer 
initial data satisfying (|1.5p , then the integrality of the subsequent terms and 
the divisibility property (II. 6|) can be proved by induction using the relations 
(|2.5p . For other special properties of integer EDS we refer the reader to 
Ward's papers [33l |3l] and to the more recent works [SJ El [TJ [25l [28l l3T] . 
If we are working with an elliptic curve E over C, then the terms Wn = 

2 

a{nn)/a{n)^ of a generalized EDS are essentially just values of the the 
division polynomials of E (see Exercises 24 and 33 in chapter 20 of [35], 
and chapter II of ^7J), possibly up to a prefactor of the form 7^""' (see 
Lemma 7 in [28j). Each Wn is an elliptic function of k that can be written 
as a polynomial in the variables {\, jj) = (p(k),p'(k)) with coefficients in 
Q[52,53] (cf. Exercise 3.7 in [26], or chapter II in [E]). For our purposes, 
in order to make use of the EDS that is the "companion" of a Somos 4 
sequence, we will need to write the Wn as polynomials in a slightly different 
set of variables (Theorem 2.2 below). 

Given any Somos 4 sequence satisfying a bilinear recurrence (II. ip with 
coefficients a, /3, there is a natural companion EDS associated with it, which 
can be defined in two different ways [211 iHj . 

Definition 1. Algebraic definition of companion EDS: For a Somos 4 
sequence, compute the translation invariant according to lll.9\) . and then the 
companion EDS {Wn}n€Z is the (generalized) EDS that satisfies the same 
fourth- order recurrence (QHP with initial values specified to he 

Wi = l, W2 = -V^, W3 = -P, Wa=IV^, (2.6) 

where 

l = a^ + (5T. (2.7) 

(Note that the coefficients of the Somos 4 sequence are then given in terms 
of the first three terms of its companion EDS by jj.^l ). j 

Definition 2. Analytic definition of companion EDS: For a Somos 4 
sequence, according to Theorem 1.2 the terms An are given explicitly by the 
formula for suitable parameters g2, g3,a,b, zq, k, and then the terms 

Wn of the companion EDS are given in terms of the Weierstrass sigma 
function by the formula 7| ) with the same parameters g2,g3, i^- 

In his original memoir [33], Morgan Ward showed that an EDS admits 
two equivalent definitions, one algebraic and the other analytic. By con- 
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struction, the translation invariant for the Somos 4 sequence and that for 
its companion EDS have the same value T = p"{k), and, upon comparing 
()1.4p with (12. 4p , the algebraic and analytic definitions of the companion EDS 
are easily seen to be equivalent over C; see also the Remarks after equations 
(1.18) and (2.14) in [T3]. The introduction of the companion EDS is very 
natural in the light of the following result due to one of us with van der 
Poorten 



Theorem 2.1. The terms An of a Somos 4 sequence satisfy the Hankel 
determinant formulae 



and 



W1W2A. 



n+m+lAn—m 



WmAn+1 Wm-lAn 
Wrn+lAn WmAn-l 



Wm+lAn+2 Wm-lAn 



(2.^ 



(2.9) 



for all m,n £ 7^, where Wm o-Te the terms of the companion EDS. 



Proof. For a purely algebraic proof of (j2.8p and (|2.9p see [2T]. For an 
analytic proof based on the three-term equation for the sigma function see 
Corollary 1.2 and Corollary 1.3 in [Hj. Note that although Wi = 1 we 
have included the Wi terms on the left hand sides above to illustrate the 
homogeneity of these expressions. □ 

It is clear that a (generalized) EDS is its own companion, and Ward's 
formula (j2.5p is the special case of (|2.8p when An = Wn- Theorem 2.1 is our 
main tool for proving integrality properties of Somos 4 sequences. To begin 
with we can use it to derive a property of EDS which is apparently new. 

Theorem 2.2. Polynomial representation for elliptic divisibility se- 
quences: For a (generalized) EDS defined by the fourth- order recurrence 
( fj.3|] with initial data Wi = 1, W2 = —\foL, W3 = — /?, W4 = T^/oL, the 
terms in the sequence satisfy W2n-\ G Z[a^,/3,T] and W2n £ \/« '^[ct'^ , (3,1] 
for all n G Z. 



Proof. Since W-n = —Wn we need only consider n > 0, and then the 
proof is by induction on n. Clearly the result is true for the initial data 
Wi, W2, W3, W4, and the next terms are 

W5 = -a^I + /3^ We = - V^(J2 + a^I - 

so the appearance of is evident. Now setting n = m + 1 and n = m + 2 
in Ward's formula ()2.5p . and putting in the initial values Wi and W2, yields 

W2m+1 = {W„,fW„,+2 - iW„,+ifWm-l (2.10) 
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and 



vV2m+2 = 1= l^--LJ^j 

respectively, which correspond to weU known identities for division polyno- 
mials - see e.g. Exercise 3.7 in [26]. Each term on the right hand side of 
()2.10p is a product of either four even index terms or four odd index terms, 
and hence by the inductive hypothesis, for m > 3, W2m+i £ ^[o^^i /Sj^T]. 
Similarly for W2m+2, both terms in the numerator on the right hand side of 
(j2.1ip consist of a product of two odd index and two even index terms, and 
so lie in aZ[a^, by the inductive hypothesis; dividing out by ^/a in the 
denominator gives the required result. □ 

Remarks. Taking the Weierstrass model 

Y'^ = X^ + AX + B, (2.12) 

which is equivalent to (I1.12|) via x = X, y = 2Y, 52 = —4^, gs = —AB, 
it is known that the division polynomials corresponding to the multiples 
[n] {X, Y) are elements of Z[A,B, X, Y^] for n odd and of 2Y Z[A,B, X, Y^] 
for n even. Hence using the equation ()2.12p to eliminate Y^, and dividing 
by 2Y where necessary, gives polynomials in 7j[A,B,X] (see Exercise 3.7 
in [26) . or [T7], chapter II). Thus we see that the result of Theorem 2.2 
corresponds to a different choice of basis for the division polynomials: after 
suitable rescaling by 2Y for n even, they are polynomials in the variables 

0-2 = i6y4^ p = A^ - GX'^A - 12XB - 3X^, 

I = 2A^ + WX^A^ - {lOX^ - 8XB)A - 2X^ - mX^B + 165^ 

with integer coefficients. As far as we are aware, this is a new result about 
division polynomials. The quantities I3,T and (via p.l2p ) are themselves 
elements of B, X]. 

Theorems 2.1 and 2.2 together imply a stronger version of the Laurent 
property for Somos 4 sequences, which lead us to establish a fairly weak set 
of criteria for integrality in the next section. 

3 Proof of the strong Laurent property 

With the above results at hand, we are now able to state our main result 
concerning a strong version of the Laurent property for Somos 4 sequences. 

Theorem 3.1. The strong Laurent property for Somos 4 sequences: 

Consider a Somos 4 sequence defined by a fourth-order recurrence with 
coefficients a, (3, initial data 741,742,^3,^44, and translation invariant T as 
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in 1^1. 13\) . The subsequent terms of the sequence are elements of the ring 
TZ of polynomials generated by these coefficients and initial data, as well as 
A^^ and the quantity X = + /5T. In other words, 

Anen = Z[a,/3,I,Af\A2,A3,A4 

for all n > 1. 

Proof. Setting n = m + 1 in each of the Hankel determinant formulae (j2.8p 
and (12. 9p and substituting Wi = 1, W2 = —\/a leads to 

{Wra)^AmAm+2 - VFm+l W^m-l^^+l iN 
A2m+l — Z (•J-J-j 

and 

. Wm+2Wm-lAm+lAm+2 — WmWm+lAmAm+3 n\ 

A2m+2 - 1=^ • (O-^j 

^/aAl 

The proof that An G TZ for n > 1 then proceeds by induction, using the 
identities and (|3.2p for odd/even n respectively. For odd n = 2m+l it is 
clear from Theorem 2.2 that, on the right hand side of ()3.ip each of the terms 
and Wm+iWm-i G and by the inductive hypothesis the other 

terms in (13. ip are in TZ. Similarly, for even n = 2m+2, on the right hand side 
of (|3.2p it is clear that both Wm+2Wm~i/\^ and WmWm+i/\^ G ^[a, P,T], 
and the result follows. □ 

Remarks. From (11131) it is clear that T e Z[a, P, Af^ , Af^ , Af^ , A^'^ ], 
so I e Z[a,P,Af^,Af^,Af^,Af^]. Similarly, applying Theorem 3.1 to 
the reversed sequence A* = ^5_n starting from A* = A^ gives An £ 
Z[a, (3,1, Ai, A2, A3, Af^] for all n < 4. Hence it follows that for all n G Z, 
the terms An G Z[a, (3, Af^ , Af'^ , A^^ , A^^], so Theorem 1.1 is a conse- 
quence of Theorem 3.1. However, note that as it stands, Theorem 3.1 is 
unidirectional (it only applies for n > 1), whereas Theorem 1.1 applies to 
both positive and negative n. 

Corollary 3.2. Integrality criteria for Somos 4 sequences: Suppose 
that a Somos 4 sequence is defined by the recurrence ( fi.i)) for coefficients 
a, (3 G 7j, and initial values Ai = ±1 with non-zero A2, A^, G Z. If [3T is 
an integer, where T is the translation invariant given by 1^1. 13\) . then An £ 7L 
for alln > 1. If it further holds that the backwards iterates Aq, A^i, A_2 G 
then An G Z for all n G Z. 

Proof of Corollary 13.21 For n > 1 the integrality is obvious, since 2 = 
o? + /3T and a, /3 G Z. If the three backward iterates adjacent to A\ are 
also integers, then to get integrality for n < 1 it suffices to apply Theorem 
3.1 to the reversed sequence Bn = A2-n, and then Bn G Z for n > 1 and 
the result follows. □ 
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Corollary 3.3. Given a,(3 G Z and T e Q with f3T G 7L, suppose that the 
quartic equation 

+ a {su^ + t^v) + P t^u^ = T stuv (3.3) 

has a solution of the form {s,t,u,v) = {Ai, A2, A^, A^), with A\ = ±1 and 
non-zero A2^ A^^ Ai^ G %. Then provided that the orbit of this set of initial 
data under is non-periodic, it produces infinitely many integer solutions 
of the Diophantine equation liS.cl]) . 

Remarks. Note that the equation ()3.3p is homogeneous in {s,t,u,v), 
so it is really the solutions with gcd{s,t,u,v) = 1 that are the interest- 
ing ones. For the sequence ()1.2p it is easy to show by induction that 
gcd(j4„, An+i, An+2, An+3) = 1 for all n; but in general for sequences where 
this holds for n = 1 then it need not be so for all n, particularly when 
gcd(a, /3) 7^ 1; see Example 13.41 for instance. Nevertheless, in such cases 
it can still be often be checked that one does get infinitely many distinct 
solutions satisfying gcd{s,t,u,v) = 1. 

Proof of Corollary 13.31 From Theorem 1 1.2 1 we see that the equation (|3.3p 
is just a rewriting of (II. 9p . Since T is a conserved quantity for (jl.ip . for 
all n the quadruple {s,t,u,v) = {An, An+i, An+2, An+3) lies on the quartic 
threefold defined by ()3.3p . and by Corollary l3.2l this is a quadruple of integers 
for n > 1. As long as the orbit of the initial quadruple (^1,^2,^3,^4) is 
not periodic (which would correspond to P G E being a torsion point), these 
quadruples are all distinct. □ 

Example 3.4. Consider the sequence defined by 

An+iAn = ll3^„+3^„+l + 90 • 113^2^2 

with initial values Ai = 1,A2 = 3, A3 = 11^,^4 = 11^ • 7 • 19; this is just 
the example previously mentioned with 1^1. 8\) . This sequence extends in both 
directions, thus: 

2498287, 1221, 7, 1, 3, 121, 177023, 2460698229, . . . (3.4) 

In this case we have integer coefficients a = 1331, /? = 119790, and with 
n = 2 in the formula lll.y\) from Theorem we calculate T = 869, so 
certainly f3T E Z. Also we have Ai = 1, and from ^3.41 ) we see that on 
each side of this, the three adjacent iterates (A2, A3, A4^ and ^0,^-1,^-2 
respectively) are all integers, so all the conditions of Corollary \3.2\ are met 
and this sequence consists entirely of integers. It is easy to see that the terms 
of the sequence grow rapidly, such that log^,i ~ Cn^ with C ~ 1.5. So the 
sequence cannot be periodic, and hence (by Corollaru \3. 3)) it yields infinitely 
many solutions of the Diophantine equation 

s^v^ + 1331(sn3 + t^v) + 119790t^u^ = m^stuv. 
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The preceding example admits a rather broad generahzation, as follows. 

Example 3.5. Pick three integers a,d,e, and choose an ordered pair of 
integer factors b, c such that 

a^d + e^ = bc. (3.5) 

Now take the Somos 4 recurrence of the form 

An+iAn = e^ An+sAn+i + adc^ {An+2f, (3.6) 

with the initial conditions Ai = 1, A2 = a, A3 = e'^ , A4 = ce^. This extends 
in both directions as 

e^ib'^d + e{b + d)), ae{b + d),b, 1, a, e^, ce^, ae^(c + de),... (3.7) 

so that the three adjacent iterates both to the left and to the right of A\ are 
integers. Thus we have a = &" , /3 = ade^ G Z and from \1.9^ we calculate 
f3T = de^{a^ + be + c) £ Z, which means that Corollary \3.S\ applies and 
An £ TL for all n£TL. 

Remarks. The preceding example can be interpreted in several different 
ways. Working over C, the equation (|3.5|) defines an affine variety V in C^, 
and then by Theorem 11.21 there is a family of elliptic curves fibred over V , 
with the sequence (j3.7p corresponding to a sequence of points on each fibre. 
Note also that by setting X = —ad, Y = de, D = bcd^ we get 

= X^ + D, (3.8) 

so that the variety V itself admits a fibration with all generic fibres {bed 7^ 0) 
being isomorphic over C to the same elliptic curve (since all the curves (13. 8p 
are isomorphic for Z) 7^ 0). Alternatively, we can consider the function field 
generated by the variables a,b,c,d,e subject to the relation (13.50 . In that 
case, Theorem 13.11 implies that the recurrence (|3.6|) generates a sequence 
of polynomials in the ring Z[a, 6, c, d, e]/ ~, where ~ is the corresponding 
equivalence relation. Note also that Example l3.4l is the particular case a = 3, 
b = 7, c = 133, d = 30, e = 11, while the original Somos (4) sequence given 
by (jl.2p is the case a = 1, b = 2, c = 1, d = 1, e = 1. Generically (provided 
that b ^ c and e 7^ 1) a different ordering of the factors b, c in (13. 5p not only 
gives a different sequence ()3.7p but also corresponds to a sequence of points 
on a different (non-isomorphic) elliptic curve. 

We should point out that the sufficient criteria for integrality in Corollary 
13.21 are not necessary conditions. To see this, it is enough to consider the 
fact that Somos 4 sequences are invariant under the two-parameter abelian 
group of gauge transformations defined by 

An — >An = db'^An, d,b £ C* , (3.9) 
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in the sense that An satisfies the same bihnear recurrence (ll.ip as does An- 
Furthermore applying (j3.9p leaves the the translation invariant T the same, 
and hence the curve E in Theorem 11.21 is preserved by the action of this 
gauge group. Moreover, if A„ E Z for all n and the transformation (j3.9p is 
applied with a, 6 G Z*, then integrality is preserved for n > 1, and if we fix 
6 = 1 then An £ Z for all n, and \An\ > 2 for \d\ > 2. 

As well as the gauge transformations (|3.9p . Somos 4 sequences exhibit 
an additional covariance property: under the action of the group of trans- 
formations given by 

An^An = ff''An, C £ C* , (3.10) 

the terms An of the transformed sequence satisfy a Somos 4 recurrence, 
which is of the same form (jl.ip but with the coefficients and translation 
invariant rescaled according to 

a = c^a, f3 = c^f3, f = c^T, (3.11) 

so that the j-invariant (jl.lOp is preserved. For each c, the transformation 
(j3.10p just corresponds to homothety acting on the Jacobian: z E Jac(£') is 
mapped to z = c~^z E Jac(S). It is clear that applying the transformation 
(|3.10p with c E Z* preserves integrality of Somos 4 sequences, as does the 
transformation 

A > /it — ( A, \^^~'^^'^~^ A 

which inserts a 1 at the A:th term of the sequence. All transformations of 
this type, that include elements of the group (I3.10p with c E Z*, have the 
effect of increasing the discriminant of the associated curve E when |c| > 1. 

Corollary 13.21 only applies when one of the terms in a given sequence 
is ±1 and sufficiently many adjacent terms are integers. If this is not the 
case, then one might try to engineer it to be so by applying a combination 
of the transformations (j3.9p and (j3.10p . in order to minimize the size of the 
discriminant as much as possible while preserving the requirement that a, (5, 
T E Z as well the the integrality of a given set of adjacent terms. However, 
without the requirement that one of the terms of the sequence should be 
±1, we can present a different set of sufficient criteria for integrality. 

Theorem 3.6. // a, /3 E Z and eight successive terms A^2i ^-i ■ ■ ■ , ^4; ^5 

are integers, and if 

gcd(a,/3) = gcd(Ai,^2) = gcd(a, ^o, -^2) = gcd(a, ^1,^3) = 1, 

then the terms An of a Somos 4 sequence are integers for all n G 

Proof. Suppose that a prime p appears in the denominator of some term, 
and suppose that p IA\. Then by Theorem 3.1, p divides the denominator 
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of T, and hence it must divide the denominator of T G Q. Setting n = 1 
and m = 3 in the result (j2.8p from Theorem 2.1, we have 



by the integrahty assumptions, whence p\a] and then p //3 because gcd(a, (5) = 
1. Thus, making use of the recurrence (II. ip gives 

A.iAs = aAoA2 + P{Aif = P{Aif ^ mod p 

But from the expression (jl.l6p with n = and n = 2 we have that p|j4_ij4oj4i 
and respectively. It follows that p divides gcd(^0;^2) and hence 

gcd(Q, ^2)- But gcd(a, Aq, ^2) = 1, contradicting the initial assumption 
that p J(Ai. 

Now by repeating the above argument with each An replaced by An+i, 
it is clear that if p is a prime which appears in the denominator of some 
term, and if p is coprime to A2, then p divides gcd(Q;, ^1, ^3). But then, 
since gcd(Q, Ai, A3) = 1, we must have also that p\A2- Thus we require that 
p\ gcd(>li,^2) and the result follows by contradiction. □ 

It has been known for some time that Somos 5 sequences, defined by 
bilinear recurrences of the form 

Tn+3''"n-2 = 0!Tn+2Tn-l + /? r„+iT„, (3.12) 

are also related to sequences of points on elliptic curves [371 UHl EI] (see 
also the unpublished results of Elkies quoted in [2j ) . The complete solution 
of the initial value problem for the general Somos 5 sequence, in terms of 
the Weierstrass sigma function, was first presented in It turns out 

that Somos 5 sequences have two algebraically independent analogues of 
the translation invariant, but for our purposes here we will only need the 
quantity J from [T3] given in terms of five initial data by 

T3''"2 \T4T1 T5T2J TsTi 

Rather than summarizing all the results of [H], we will just state our main 
theorem concerning the Laurent property and refer to some of these results 
as needed in the proof. 

Theorem 3.7. The strong Laurent property for Somos 5 sequences: 

For a Somos 5 sequence defined by a fifth-order recurrence 1^3. 12\) with coef- 
ficients a, (3 and initial data ri, T2, r3, r4, rs, the quantity 

J ^Tn+2rn-l ^^f TnTn-l ^ K)' (3.14) 

TnTn-l Tn+lTfi \'''n+lTn-2 Tn+2Tn-l J ''"n+2Tn-2 
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is independent of n. The subsequent terms of the sequence are elements of 
the ring of polynomials generated by these coefficients and initial data, as 
well as , and the quantity X = (3 + aJ. In other words, 

for all n> 1. 

Proof. The fact that J given by (j3.14p is independent of n is part of The- 
orem 2.5 in [H], where it is rewritten in terms of the quantity 

According to Corohary 2.10 in [14J (or the comments in Section 7 of |21j). 
each Somos 5 sequence also has a companion EDS associated with it, with 
terms denoted a„ there, such that the analogue of (|2.9p . given by 



(3.15) 



holds for all m, n G Z. Using further results of that Theorem 2.5, we can 
define this companion EDS by the Somos 4 recurrence 

~2 ~ t \2 

with initial data ai = 1, 02 = —jl, 03 = d, 04 = where 

ii = + ajf'^. (3.16) 

It follows from Theorem 2.2, making the replacements y/a jl, (3 ^ —a, 
T — > /3, that a2n+i and a2n/A are both elements of Z[/i^, d, /?], so by (I3.16p . 
noting that T = /i^ we have a2n+i G ^[dj/^jX] and a2n G /^^[o^j for all 
n. Setting n = m + 1 and n = m + 2 in (|3.15p yields two identities for r„ 
which for even/odd index n have respectively ri and T2 in the denominator. 
Thereafter the proof proceeds similarly to that of Theorem 13.11 □ 



Remarks. Similarly to the Remarks after Corollary 13.21 the Laurent prop- 
erty for Somos 5 sequences, as proved in [8], is a consequence of Theorem 
13. 7i Also, the latter has two immediate corollaries which are the respective 
analogues of 13.21 and 13.31 

Corollary 3.8. Integrality criteria for Somos 5 sequences: Suppose 
that a Somos 5 sequence is defined by the recurrence I13.12\) for coefficients 
a,f3 £ 7j, and initial values ti = ±1, T2 = ±1, with non-zero t^jT/^jT^ G Z. 
If the quantity aJ is an integer, with J given by 1^3. 13\) . then Tn £ 1' for all 
n > 1. If it further holds that the backward iterates ro,T_i,r_2 € Z, then 
Tn € for all n G Z. 
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Corollary 3.9. Given a, /3 G Z and J G Q with aJ G Z, suppose that the 
quintic equation 

{sw + a u^){sv'^ + Pw) + (3 tu^v = J stuvw (3-17) 

has an integer solution of the form {s,t,u,v,w) = (ti, r2, rs, r4, rs), with 
\ti\ = 1, \t2\ = 1 and non-zero t^^t^^t^ E Then provided that the orbit 
of this set of initial data under 1^3. is non-periodic, it produces infinitely 
many integer solutions of the Diophantine equation ^3.17^ . 

Example 3.10. Take a = 14641 = 11^, /3 = 1771561 = 11^ in (TTM> with 
initial data 847, 8, 1, 1, 33. From ^3.14^ we calculate J = 627, hence aJ S Z, 
and the sequence extends on either side of these values as 

805255, 847, 8, 1, 1, 33, 6655, 19487171, . . . , (3.18) 

so all the conditions of Corollary \3.8\ hold. Thus Tn £ Z for all n £ 1^. By 
Corollaru \3.y\. for these values of a, (3, J the sequence i3.18\) gives infinitely 
many quintuples of integer solutions of the corresponding Diophantine equa- 
tion, of the form 1^3.17 ). 

Remarks. It is shown in [Tl] (Proposition 2.8) that the subsequences of 
a Somos 5 sequence consisting of the even/odd index terms respectively 
satisfy the same Somos 4 recurrence relation. If we index the terms of the 
preceding example such that ti = 8, T2 = 1, T3 = 1, T4 = 33 etc. then 
we find that the even index subsequence of ()3.18p is essentially the same 
as the sequence (13. 4p in Example 13.41 up to applying transformations of 
the form (|3.9|) and (|3.1U|) . To be precise, we have T2n = 11^^"^ An where 
Q(n) = (n — l)(3n — 4)/2. Indeed, both sequences correspond to a sequence 
of points on the same elliptic curve with j'-invariant (from p.lOp ) 

5^ • 236 • 1013^ 
~ 24.117.172.37.1069' 

It is clear from the denominator that the curve has bad reduction modll, 
and all terms of the sequence reduce to mod 11 apart from the central 
terms 8, 1, 1. 

Note that so far all of our criteria for integrality have been stated for 
Somos recurrences with integer coefficients. However, even when a, /3 G Z 
the associated curve (|1.12p must be defined over Q, or more generally over 
Q{Vd() where a is the square-free part of a, in order to consider the sequence 
of points (|2.3p . Hence it is natural to define Somos sequences over Q, which 
leads to the question of whether integer sequences can arise for non-integer 
a OT f5 € Q. The following example illustrates that this is indeed the case. 

Example 3.11. Consider the Somos 4 sequence defined by ( fi. j|) with a = 
— 1/2, (3 = 1 and Ai = 1, = —2, = 2, A3 = 1. This sequence has terms 

. . . , 2, 10, -4, 2, 3, 1, -2, 2, 1, 5, 2, 12, -26, 34, 236, 352, -1912, . . . , (3.19) 
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and consists entirely of integers. 



The sequence (I3.19P was presented in the Robbins (bilinear) forum [22] 
by Michael Somos, who found that it should correspond to a sequence of 
points on an elliptic curve E with j-invariant 

33 • 190513 

this was confirmed by Elkies 0- Although Theorem 13.11 does not apply di- 
rectly in this case, we show in the next section how a slight modification 
of the methods used previously provides a proof of integrality for a one- 
parameter family of Somos 4 sequences that generalizes Example 13.111 



4 One-parameter family of integer sequences 

In this section we consider a one-parameter family of sequences 

. . . , 1, -N, N,1,1 + N'^,N, N^ + 2N, -N^-2N'^-2, + + . . . (4.1) 

which is defined by the Somos 4 recurrence (jl.ip with coefficients and initial 
data given by 

a = -l/N, P = l, ^1 = ^4 = 1, A2 = -As = -N. (4.2) 

When N = 2 this is the sequence (I3.19P in Example 13.111 above. The in- 
tegrality properties of the sequence (|4.ip do not follow from Theorem 13. H 
because the parameter a is non-integer for integer N ^ ±1. Nevertheless, 
we are able to prove that the sequence (14. ip consists entirely of polynomials 
in N with integer coefficients, and so the integrality of this sequence when 
N £ Z follows immediately. Prom ()1.9p we compute the translation invari- 
ant T = —N'^ — 1/A^^, and hence from (|2.7p we have I = —N'^. Thus we see 
that X E Z for N £ Z, which leads to special properties for the companion 
EDS associated with (fiT]) . 

Theorem 4.1. For the companion EDS associated with j^^. j[ ), defined by 
the fourth order recurrence 1^1. 3\) with initial data Wi = 1, W2 = —iN~^/'^, 
W3 = —1, VF4 = —iN^/"^, the terms in the sequence satisfy W2n-i £ Z[A^^], 
Win e iN^^'^ Z[iV4] and W4„+2 e «A^"^/^ I'[N'^] for allneZ. 

Proof. This follows the same pattern as the proof of Theorem 2.2, using 
the two identities (pTO]) and ([2TT|) . except that in I^JT^ it is necessary to 
consider the cases of odd/even m separately. □ 

^See|http : //www. math. wise . edu/~propp/somos/elliptic| 
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Theorem 4.2. Consider the Somas 4 sequence defined by the recurrence 
U.l\) with coefficients a = —1/N, /3 = 1 and initial data Ai = 1, A2 = —N, 
^3 = A^, ^4 = 1. This sequence satisfies An E Z[A^] for alln € Z. Moreover, 
An £ 'Z[N'^] whenever n = or 1 mod 4, while An S NZ,[N'^] whenever 
n = 2 or 3 mod 4. 

Proof. This proceeds by induction for n > 1, analogously to the proof of 
Theorem 3.1, using the two identities ()3.ip and (13. 2p . except that each of 
these formulae must be considered separately for odd/even m to get the 
different properties of An as n varies mod 4. We omit further details, except 
to mention that to cover n < 1 it suffices to consider the reversed sequence 
An = A^^n and then proceed as before. □ 



Corollary 4.3. For all G Z, the sequence j[ j defined by the Somos 4 
recurrence with coefficients and initial data given by ^4-^ is an integer 
sequence. 

Thus we see that the sequence p.l9p is the particular case = 2 of the 
family of integer sequences given by (14. ip with N ^TL. It is then interesting 
to consider the family of elliptic curves corresponding to these sequences. 
By Theorem 1.2, using the formulae ()2.1|) . (12. 2p and p.3p from section 2, 
for each 7^ the sequence (|4.ip is associated with the sequence of points 



(iV^ - 1)2 A^ 2iN^ \ ( {N^-\f i 

+ \n\ 



\ 12Ar3 (Ar2 - 1)2 ' (iV2 - 1)3 i ' '\ 12iV3 ' ViV 
on the curve 



with invariants 



53 (iV) = 
and j-invariant 



52 (iV) 



y2 = 4x3-52(A^)x-g3(iV) (4.3) 



12A^6 

^24 _ gjY20 ^ 51^16 _ 5Q^i2 _^ ^95jy8 + 30Ar4 _^ 
216iV9 



[N^^ - 4A^i2 ^ 30^8 ^ 20Af^ + 1)^ 
' ~ Ari6(Ari2_57v8 + 39Ar4 + 2) ' ^ ' 

(The latter reproduces the value (j3.20p when A^ = 2.) 

Strictly speaking the recurrence for ()4.ip does not make sense when 
A^ = 0, but the sequence (|4.ip can still be defined by setting A^ = in each 
of these polynomials. It can be seen that when A^ = the values of these 
polynomials are given by 

A^m = (-1)™~12™(— A^n.+1 = 2-(— ^4„_i = A^m+2 = 

(4.5) 
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for all m. The case = 1 is also special, because in that case the base point 
Q = oo, and we have the sequence 

1, -1, 1, 1, 2, 1, 3, -5, 7, 4, 23, 29, 59, -129, 314, . . . (4.6) 

with ^0 = 0, which corresponds to the multiples [n]P of the point P = (0, i) 
on the curve 

= Ax^ - Ax - 1, (4.7) 

for which j = 2^^ ■ 3^/37. In that case, the sequence (14. ID is almost identical 
to its companion EDS; more precisely, we find that An = {—i)^~^Wn for all 
n. Over C (sending x — > —x, y ^ iy), the curve (14. 7p is isomorphic to the 
curve associated with the Somos (4) sequence ()1.2p as found in [13j, which 
is not surprising when one observes that the terms of ()1.2p are precisely the 
odd index terms of the sequence (j4.6p . 

The equation (|4.3p can be thought of as defining a curve E/C{N). Let- 
ting vary this gives an elliptic surface fibred over P^: with the exception 
of a finite set of values, each value of defines a non-singular elliptic curve 
(see Chapter III in [27]). The singular fibres correspond to the zeros of the 
denominator of j{N), i.e. N = and the zeros of iV^^ - 5N^ + 397V^ + 2, as 
well as the fibre at infinity. If we set S = N^, then it is clear from (14. 4|) that 
j{N) € C(5). It turns out that over C(\/iV), the rational elliptic surface 
()4.3p is birationally equivalent to the Weierstrass model 

f + 8Sy = x'^ + {S - if x^ - 8S{S + l)x (4.8) 

which is defined over C{S). In terms of the model ()4.8p the sequence (14. ip 
corresponds to the sequence of points 

Q+[n]P= ( 8SVS V^^j/ \ 

^ ^ ^ (V5-i)2 {Vs-if J ^ ' ; 

defined over the extended function field C(\/S')- 

5 Conclusions 

Using the connection between Somos 4 sequences and their associated com- 
panion EDS, together with the identities in Theorem 2.1, we have shown 
that these sequences possess a stronger variant of the Laurent property in 
[S] , which is expressed by means of the conserved quantity Z, given by the 
formula ()2.7p in terms of T, the translation invariant ()1.9p . This strong 
Laurent property has produced a set of sufficient criteria for integrality of 
Somos 4 sequences with integer coefficients, and we have given a similar set 
of criteria for Somos 5 sequences as well. 
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The fact that every Somos 4 sequence admits a conserved quantity rests 
on the hnk with shifts by muhiples [n]P of a point P on an elhptic curve E 
|31j . or equivalently with a discrete Unear flow on the Jacobian Jac(S). This 
can also be understood from the connection with integrable maps (32] : via 
the subsitution /„ = An+i An-i / {An)'^ , the bilinear recurrence (|l.ip yields 
the second-order difference equation 

fn+l{fnffn^l = afn + P, 

which has the conserved quantity (jl.l4p . preserves the symplectic form uj = 
(/n-i/n)~^ dfn-i^dfn, and is a degenerate case of the family of maps studied 
by Quispel, Roberts and Thompson [23j. 

There is another aspect of the Laurent property, that was stated as a con- 
jecture in which is that for the octahedron recurrence (also known as the 
discrete Hirota equation [36], an integrable partial difference equation) the 
corresponding Laurent polynomials in the initial data and coefficients have 
all positive coefficients, and furthermore all these coefficients are 1. Speyer 
has given a combinatorial proof of this conjecture, using perfect matchings of 
certain graphs [30j, and Speyer and Carroll have also used combinatorics to 
prove the analogous result for the four-term cube recurrence (also known as 
the Hirota- Miwa equation). Since the Somos 4 recurrence can be obtained as 
a reduction of the discrete Hirota equation, Speyer notes that his result im- 
plies that the Laurent polynomials in Z[a, /3, ^J^] for Somos 
4 have all positive coefficients. Whether this result has any interpretation 
in terms of elliptic curves is not clear. 

Although our most general integrality criteria in Section 3 were for- 
mulated for Somos 4 recurrences with integer coefficients, the family of se- 
quences (|4.1|) illustrates the fact that a and (3 need not be integers for integer 
sequences to result. The sequence of polynomials ()4.ip provides further in- 
sights into the integer sequences corresponding to G Z. In particular, 
consider the case of integer N > 2, and suppose p > 2 is a prime with p\N. 
In that case, due to the factor of A^^^ in the denominator of the j-invariant, 
the curve ()4.3p - or equivalently (14. 8|) - has bad reduction at the prime p. 
Furthermore, by Theorem 4.2 it is clear that A4m+2 = ^4m+3 = mod p, 
while from the values ()4.5p of these polynomials at A^ = it follows that 
P A^im and p /\A^m+\ for all m. This means that the distance, or gap 
length [24J, between successive multiples of the prime p appearing in the 
sequence alternates between one and three. This is in contrast with Robin- 
son's numerical observations for the Somos (4) sequence (|1.2p . which led him 
to the conjecture that for this particular sequence, the gap length should 
be constant for any prime p |24j . The p-adic properties of EDS have been 
considered recently by Silverman [28], and we expect that similar consid- 
erations will lead to a better understanding of prime divisors in Somos 4 
sequences. 
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Certain higher order Somos sequences have also appeared in connection 
with division polynomials [3], higher order integrable maps [1] and divisor 
sequences for hyperelliptic curves [iHlEO]. Based on naive counting argu- 
ments, in [13] one of us made the conjecture that all Somos sequences should 
arise from appropriate divisor sequences on such curves. In the forum |22j . 
Elkies has given much stronger arguments to the contrary, based on a con- 
jectured theta function formula for the iterates. Moreover, the Laurent 
property is failed for Somos 8 recurrences; and for the sequence generated 
by the recurrence 

with initial values Sj = 1 for j = 1, ... ,8 the logarithmic heights of the 
rational iterates have exponential growth, so that log h{Sn) ~ Kn as n ^ oo 
with K « 0.23; this growth is incompatible with a discrete linear flow on an 
Abelian variety. Recently, Halburd has proposed that polynomial growth of 
logarithmic heights should be an integrability criterion for birational maps 

One of us has found that for Somos 6 (and also Somos 7) recurrences, 
there are two independent conserved quantities, analogous to the translation 
invariant ()1.9p for Somos 4. For these recurrences of sixth and seventh 
order, it is possible to given an explicit formula for the iterates in terms of 
two-variable theta functions, or equivalently in terms of the Kleinian sigma 
functions for an associated curve of genus two. Recently, Kanayama has 
derived multiplication formulae for genus two sigma functions [15]. One of 
us has pointed out that the statement of Proposition 3 in |15j is incorrect, 
but Kanayama has shown us a corrected version of this result [TH] , which we 
have used to derive the solutions of Somos 6 and Somos 7 recurrences. The 
full description of the genus two case is the subject of further investigation. 
Acknowledgements. AH is grateful to James Propp, Noam Elkies and 
other members of the Robbins bilinear forum [22] for helpful correspondence 
on related matters. 
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